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Disclaimers 
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INTRODUCTION 


The  quantitative  measurement  of  the  absorption  coefficient  of 
atmospheric  dust  is  a  problem  made  difficult  by  the  complicated  nature 
of  the  material.  Since  dust  samp I es  consist  of  mixtures  of  finely 
powdered,  highly  scattering  materials,  conventional  transmission 
spectroscopy  is  not  directly  applicable.  The  basic  problem  is  to 
determine  what  fraction  of  the  incident  light  attenuated  by  a  sample  of 
dust  is  lost  due  to  scattering,  and  how  much  is  absorbed  by  the  particles 
themselves.  The  measurement  of  absorbed  energy  must  then  be  related 
to  the  basic  optical  absorption  coefficient  of  the  sample.  This  latter 
problem  is  not  trivial  since  the  Bouguer-Lambert  law,  which  normally 
solves  this  problem  for  nonscattering  media,  is  not  directly  applicable 
here. 


In  spite  of  the  difficulties  in  making  such  measurements,  several 
workers  have  devised  methods  for  solving  the  problem,  with  varying  degrees 
of  success.  Fisher  [1 1 D  has  obtained  resuits  based  on  a  diffuse  transmittance 
technique  and  the  assumption  that  the  Bouguer-Lambert  law  can  be  used. 

Volz  [2,  3U  has  appl ied  variations  of  the  potassium  bromide  pressed  disk 
method,  commonly  used  for  qualitative  work  in  the  infrared  region,  to 
gain  some  information  about  the  absorption  coefficient  of  dust.  Grams 
et  a  I .  (X)  have  estimated  the  imaginary  refractive  index  of  airborne  fly 
ash  by  inferring  the  value  needed  to  account  for  their  measured  values 
of  laser  radar  backscatter  cross  section  and  particle  size  distribution. 
Previous  work  in  this  laboratory  D>3  has  developed  a  somewhat  complicated 
method  for  applyinc  the  Kubelka-Munk  theory  of  diffuse  reflectance 
to  the  problem.  T! is  report  describes  a  straightforward  method  of 
obtaining  the  absoption  coefficient  from  a  single  diffuse  reflectance 
measurement.  The  ~inal  section  is  a  discussion  of  the  results  obtained 
here  in  compariso1  to  those  obtained  by  other  workers. 

THEORY 

In  this  work  we  regard  a  sample  of  atmospheric  dust  as  a  powdered  material 
of  unknown  shape,  size  distribution,  refractive  index  and  absorption 
coefficient.  The  object  is  to  infer  the  absorption  coefficient,  or 
equivalently,  t^e  imaginary  part  of  the  complex  refractive  index  of  the 
material,  from  'leasurement  of  some  optical  property  of  the  sample. 

The  most  stra i-.htforward  measurement  that  contains  information  about 
the  absorptior,  coefficient  of  a  powder  is  a  measurement  of  its  diffuse 
reflectance,  ?  .  It  is  well  known  that  the  diffuse  reflectance  of  a 
powder  is  hig>Ty  dependent  on  the  absorption  coefficient.  R^  is  high  when 
the  absorpt ic. n  is  low,  and  of  course  is  low  for  strongly  absorbing  materials. 
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For  this  reason  the  diffuse  reflectance  of  a  powdered  sample  of  a  solid 
material  plotted  as  a  function  of  wavelength  has  the  character  of  a 
transmission  spectrum  cf  the  material  itself.  The  resemblance  is 
only  qualitative,  however. 

To  obtain  qualitative  information  from  a  diffuse  reflectance  nnasure- 
ment,  it  is  necessary  to  have  some  analytical  relationship  between  the 
diffuse  reflectance  and  the  desired  optical  property  of  the  powder,  in 
this  case  the  absorption  coefficient.  Several  theories  of  diffuse  reflect¬ 
ance  provide  such  a  connection.  Some  are  very  satisfactory  quantitatively, 
but  require  considerable  a  priori  knowledge  about  the  sample,  such  as 
size  or  real  refractive  index.  For  the  work  reported  here  the  Kubel  ka-NVjnk 
theory  was  used,  since  it  requires  no  other  knowledge  about  the  sample 
at  all.  The  price  of  this  convenience  is  the  assumption  that  the  powder 
is  an  isotropically  scattering  medium,  a  point  we  will  come  back  to  later 
in  this  section,  and  a  definition  of  absorption  coefficient  slightly 
different  from  that  oommonly  used  in  the  Bouguer-Lambert  law. 

The  Kubelka-Munk  theory  as  it  is  applied  in  this  work  is  well  known  and 
will  not  be  discussed  in  detail  here.  The  reader  is  referred  to  two 
excellent  books  [6,  7 3  on  the  subject  for  details  on  its  origin  and 
derivation.  The  equations  used  here  and  the  notation,  as  much  as  possible, 
are  taken  from  the  work  by  Kortum  £7],  The  nota+ion  is  summarized  in 
Table  I. 


K 

I 


The  Kubel ka-Munk  theory  treats  a  powder  sample  as  a  continuous,  turbid 
medium.  In  a  layer  of  infinitesimal  thickness,  fractions  s  and  k  per 
unit  path  length  of  the  incident  light  are  scattered  and  absorbed, 
respectively.  The  quantities  s  and  k,  in  units  of  cm  ,  are  called  the 
Kubel ka-Murk  scattering  and  absorption  coefficients  It  is  known  that 
the  coefficient  k  is  very  closely  related  to  the  absorption  coefficient 
in  the  familiar  Bouguer-Lambert  law.  Discussions  of  this  point  are 
available  in  the  references  [6,  7,  8],  The  object  of  the  work  reported 
here  is  to  determine  the  value  of  k  for  a  sample  of  atmospheric  dust,  and 
to  infer  the  imaginary  refractive  index  from  it. 

One  c*  the  most  important  accomplishments  of  the  Kubelka-Munk  theory 
is  an  equation  relating  the  diffuse  reflectance  of  a  powder  to  the 
scattering  and  absorption  coefficients: 
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TABLE  I 

DEFINITIONS  OF  SYMBOLS 


c  Concentration,  In  males  per  liter,  of  a  pure  powder  sample 

c*  Concentration,  In  males  per  liter,  of  sample  In  a  mixture  of  sample 

and  diluting  agent 

d  Thickness,  in  cm,  of  a  sample  layer 

k  Kubelka-Munk  absorption  coefficient  in  cm"* 

k*  Kubelka-Munk  absorption  coefficient  in  cm  *  for  a  mixture  of  sample 

and  diluting  agent 

n’  The  imaginary  part  of  the  complex  refractive  index  of  atmospheric  dust 

R  Diffuse  reflectance  of  a  layer  of  powder  of  thickness  d 

Rg  Diffuse  reflectance  of  the  flat  black  paint  used  to  coat  the  bottom 

of  a  shallow  aluminum  sample  dish 

R^  Diffuse  reflectance  of  a  layer  of  powder  with  a  thickness  large 
enough  that  no  light  passes  through 

R*  Diffuse  reflectance  Rw  measured  for  a  mixture  of  sample  and  diluting 
agent 

s  Kubelka-Munk  scattering  coefficient  in  cm"* 

s*  Kubelka-Munk  scattering  coefficient  i n  cm  *  for  a  mixture  of  sample 

and  diluting  agent 

W  The  weight  of  pure  atmospheric  dust  sample  required  to  completely 
fill  the  vol jme  of  the  sample  dish 

W*  The  actual  weight  of  atmospheric  dust  in  the  sample  dish  when  R*  is 
measured 


This  expression  is  known  as  the  Kubelka-Munk  remission  function.  This 
provides  a  determination  of  the  ratio  k/s  from  one  measurement  of  refiee"- 
ance,  R  .  In  order  to  arrive  at  values  of  k  and  s  separately,  more 
information  is  required. 

The  reflectance  R  ,  by  definition,  is  the  diffuse  reflectance  of 
an  "infinitely  thTck"  layer  of  the  powder;  this  means  a  layer  so  thick 
that  no  detectable  light  is  transmitted  through  the  sample.  In 
practice  this  turns  out  to  be  a  thickness  of  a  few  millimeters,  depend !«  *; 
on  the  value  of  k  and  s  for  the  material.  Suppose  that  we  measure  Rm 
on  a  suitably  thick  layer  of  sample,  and  then  also  measure  R,  the 
diffuse  reflectance  of  a  much  thinner  layer  of  thickness  d.  This 
con  be  done  by  placing  the  powder  in  a  shallow  dish  of  dep+h  d.  Suppose 
also  that  we  have  made  the  bottom  of  the  dish  nearly  black,  and  have 
previously  measured  its  diffuse  reflectance  Rg.  Kortum  [7,  9]  has  she-  ■ 
from  the  Kubel ka-Munk  theory  that  s,  Rw,  R,  Rg,  and  d  are  related  by 
the  expression 
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and  that  s  can  be  determined  satisfactorily  from  measurement  of  these 
quantities.  Having  obtained  s  in  this  manner,  one  can  then  solve  for 
k  in  Eq.  (I).  This  provides  a  means  of  obtaining  the  absorption  coefficient 
of  a  powder  from  three  reflectance  measurements  and  one  thickness  measure¬ 
ment. 

The  method  described  above  has  been  used  successfully  to  obtain  k 
and  s  for  very  weakly  absorbing  (nearly  white)  powders.  But  it  is  not 
a  satisfactory  solution  to  the  problem  of  interest  here  for  three  reasons. 
One  is  that  for  a  sample  of  atmospheric  dust,  the  absorption  coefficient 
is  relatively  high,  and  therefore  even  a  thin  layer  appears  "infinitely 
thick,"  making  it  difficult  to  determine  R,  R^,  and  d  to  a  satisfactory 
degree  of  accuracy.  Second,  it  is  nor  easy  to  collect  a  sufficiently 
massive  sample  of  atmospheric  dust  to  carry  cut  these  kinds  of  measurements. 
The  third  problem  arises  from  the  assumption  in  the  Kubel  ka-IAjnk  theory 
that  the  sample  powder  scatters  light  isotropically.  Experiment  shows 
that  weakly  absorbing  powders  are  isotropic  scatterers  to  a  very  good 
degree  of  approximation.  But  this  is  not  necessarily  true  for  strong 
absorbers  like  atmospheric  dust. 

The  difficulties  mentioned  above  arise  from  the  fact  that  atmospheric 
dust  is  a  relatively  strong  absorber  of  light.  Therefore,  they  can  be 
avoided  by  diluting  the  dust  sample  with  some  other  powder  which  has  a 
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trivially  low  absorption  coefficient.  One  can  then  apply  Eqs.  (I)  and 
(2)  to  the  diluted  sample,  and  obtain  an  absorption  coefficient  for  the 
mixture.  The  remaining  problem  is  to  infer  the  absorption  coefficient 
of  the  pure  dust  sample  from  that  of  the  mixture  and  knowledge  of  the 
dilution  ratio. 

First  we  look  at  the  problem  of  obtaining  the  absorption  coefficient 
of  a  mixture  of  atmospheric  dust  and  some  other  powder  with  low  absorption 
coefficient  used  as  a  diluting  agent.  We  can  then  write  Eq.  (I)  in 
the  form 


(I  -  R*)r  k* 

- =_  =  —  ,  (3) 

2R  s* 

00 

where  the  asterisk  indicates  that  we  are  deal ing  with  the  properties 
of  the  diluted  sample.  Kortum  has  shown  that  If  the  sample  Is  highly 
diluted,  the  scattering  coefficient  s  of  the  mixture  is  not  affected  by 
the  presence  of  'he  sample,  so  for  a  value  of  s#  we  can  use  the  scattering 
coefficient  for  the  pure  dilutant.  This  quantity  can  easily  be  determined 
by  measuring  R^,  Ro,  Rg,  and  d,  then  applying  Eqs.  (I)  and  (2)  as 
discussed  above,  since  the  diluting  agent  is  not  a  strong  absorber  like 
the  pure  dust  sample.  We  can  measure  R£  directly  for  the  diluted  dust 
sample,  and  determine  k*  from  Eq.  (3). 

Now  we  have  a  valid  measurement  of  k*,  the  absorption  coefficient 
of  the  diluted  dust  sample.  What  we  really  want,  however,  is  the  absorp¬ 
tion  coefficient  of  the  pure  dust  sample,  undiluted  by  air  spaces  or 
other  white  diluting  agents.  Previous  work  O ,  93  has  shown  that  for 
strongly  diluted  samples,  the  absorption  coefficient  k*  is  directly 
proportional  to  the  sample  concentration,  in  moles  per  liter,  over  several 
orders  of  magnitude  of  concentration.  This  is  analogous  to  Beer's  law 
in  conventional  transmission  spectroscopy  of  solutions.  For  the  case 
of  powdered  materials,  derivations  from  linearity  occur  when  the  sample 
is  not  sufficiently  diluted,  and  the  assumption  that  the  scattering  co¬ 
efficient  of  the  mixture  is  that  of  the  pure  diluting  agent  is  not  valid. 
Nonlinearities  also  occur  for  extremely  weak  dilutions  because  the 
absorption  coefficient  of  the  mixture  is  so  small  that  the  intrinsic 
absorption  coefficient  of  the  diluting  agent  can  no  longer  be  considered 
trivial.  In  practice  for  the  kinds  of  materials  dealt  with  here,  the 
linear  region  extends  typically  over  a  dilution  range  of  one  part  sample 
in  10*  parts  diluting  agent  to  one  part  sample  in  10*  parts  diluting 
agent. 
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We  now  assume  that  we  have  measured  R*  for  a  diluted  sample  of 
concentration  which  fails  in  the  linear  region,  and  have  calculated  the 
corresponding  value  of  k*.  Then  we  can  obtain  the  absorption  coefficient 
k  for  a  pure  sample,  free  from  dilution  by  any  agent,  from  the  expression 


k  = 


(4) 


where  c*  is  the  concentration  in  moles  per  liter  of  sample  in  one  dilution, 
and  c  is  the  concentration.  In  moles  per  liter,  the  sample  would  have 
if  it  were  free  of  air  spaces  or  other  diluting  agents.  Since  concentration 
in  males  per  liter  is  directly  proportional  to  the  weight  per  unit 
volume  of  sample,  we  can  replace  the  ratio  c/c*  in  Eq.  (4)  with  the 
ratio  W/W*,  where  V#*  is  the  weight  of  atmospheric  dust  in  our  sample 
volume,  and  W  is  the  weight  of  dust  that  would  be  required  to  fill  the 
sample  volume  completely  with  the  dust.  Making  this  change  in  Eq.  (4) 
and  substituting  for  k*  from  Eq.  (3),  we  have 


s*( I  -  R*)2  W 
00 

2R*  W* 


(5) 


This  expression  gives  the  Kubelka-Munk  absorption  coefficient  for 
the  pure  atmospheric  dust  sample  as  a  function  of  one  reflectance 
measurement,  a  weight  measurement,  and  a  predetermined  scattering  coefficient 
for  the  diluting  agent. 


MEASUREMENTS 

All  of  the  diffuse  reflectance  measurements  were  made  with  a  Cary  14 
spectrophotometer  equipped  with  the  manufacturer's  25  cm  integrating 
sphere,  as  described  by  Hedelman  and  Mitchell  [103.  In  such  an  instrument 
the  quantity  measured  is  the  relative  reflectance,  that  is,  the  reflectance 
of  the  sample  divided  by  the  reflectance  of  some  material  used  as  a  reference 
standard.  Thus,  all  measurements  must  be  corrected  by  multiplying  the 
measured  reflectance  by  the  absolute  reflectance  of  the  reference  material 
to  determine  the  absolute  reflectance  of  the  sample.  In  the  work  done 
here,  highly  refined  BaSO^  [llD  was  used  as  a  reference  standard,  and  the 
absolute  reflectance  of  this  material  was  obtained  from  the  work  of  Grum 
and  Luckey  [l2j. 


8 


As  a  diluting  agent,  the  same  BaSO^  used  as  a  reference  standard  was 
chosen.  A  shallow  aluminum  dish  a  fraction  of  a  millimeter  deep,  painted 
with  flat  biack  enamel,  was  used  to  measure  R,  the  reflectance  of  a  thin 
layer  of  the  white  powder.  Rg,  the  reflectance  of  the  black  dish,  was 
also  measured-  The  scattering  coefficient  was  determined  from  Eqs.  (I) 
and  (2)  and  used  as  the  value  of  s*,  the  scattering  coefficient  of  the 
d  i  I uted  samp  I es . 

The  dust  samples  used  here  were  col  I  acted  on  the  surface  of  a  0.45  ym 
pore  size  membrane  filter  at  the  input  of  a  laboratory  vacuum  pump.  Each 
sample  was  collected  over  a  period  of  ten  days,  from  a  location  3  meters 
above  the  ground  In  a  desert  basin  In  southern  New  Mexico.  The  dust 
collected  on  the  filter  was  gently  scraped  off,  and  an  appropriate  quantity 
was  mixed  with  34  grams  of  the  BaS04  diluting  agent  Cl 3].  Mixing  of  the 
BaS04  and  dus+  was  done  'n  a  polystyrene  vial  in  a  vibrator  mixer  for  a 
period  of  5  minutes.  The  mixture  was  then  pressed  into  a  I  cm  deep  by 
5  cm  diameter  lucite  dish  using  a  clean  glass  plate.  The  reflectance,  R*, 
was  measured  directly,  and  the  weight,  W*,  of  dust  in  the  lucite  dish 
was  recorded.  The  quantity  W,  the  weight  of  dust  that  would  completely 
fill  the  dish,  was  calculated  from  the  Known  volume  of  the  dish  and  the 
assumption  that  the  average  specific  gravity  of  atmospheric  dust  materials 
is  2.4,  as  reported  by  ether  workers  lI 4U. 

Different  quantities  of  dusf  were  used  to  insure  that  the  concentration 
chosen  was  in  the  I  inear  range  as  discussed  above,  and  to  determine  over 
what  wavelength  range  the  measurement  could  be  made.  It  was  found  that 
a  suitable  linear  relationship  existed  for  wavelengths  from  0.3  to 
I. I  ym.  Outside  of  this  range,  the  intrinsic  absorption  coefficient  of 
the  BaS04  was  too  high  to  be  considered  negligible.  It  was  found  that 
quantities  of  from  3  to  30  mg  of  dust  in  the  sample  dish  produced  satisfactory 
results.  From  measurements  of  R*  at  different  wavelengths,  and  knowledge 
of  W,  W*,  and  s*  as  determined  earlier,  k  was  calculated  as  a  function 
of  wavelength  from  Eq.  (5).  Typical  results  are  listed  in  Table  II  and 
discussed  in  the  next  section. 

A  programmable  desk  calculator  was  used  to  perform  the  calculations 
required  in  this  work.  By  estimating  the  errors  in  each  of  the  input 
measurements  involved  and  determining  the  effect  of  each  on  the  computed 
value  of  k,  an  estimate  of  +he  experimental  error  was  made.  The  largest 
source  of  error  was  found  to  be  in  the  determination  of  s* ,  because  of 
the  difficulty  of  preparing  a  uniform  layer  of  powder  a  few  tenths  of  a 
millimeter  thick.  It  was  estimated  that  s*  might  be  in  error  by  as 
much  as  20%.  Photometric  errors,  assumed  to  be  about  2$,  were  found  to 
contribute  an  error  of  about  5?  to  the  determination  of  k.  Since  s*  is 
not  determined  individually  for  each  computation  ot  k,  its  uncertainty 
affects  the  absolute  accuracy  of  the  measurement,  but  not  the  repeatability 


9 


=|3sk?' 


TABLE  II 

ABSORPT I ON  COEFFICIENT  AND  IMAGINARY  REFRACTIVE 
INDEX  FOR  THREE  ATMOSPHERIC  DUST  SAMPLES 


1 

\,m 

Sample  #1 

.  -1  . 
k,  cm  n ' 

k,  cm 

Sampi e  #2 

n* 

Sample  #3 

-1 

k,  cm  n 

1 

&. 

.35 

5805 

.016 

6312 

.018 

5803 

.016 

| 

.40 

3298 

.01  1 

4239 

.013 

2837 

.009 

1  ' 

.45 

2275 

.008 

3230 

.012 

1814 

.007 

r- 

.50 

1843 

.007 

2432 

.010 

1478 

.006 

E; 

.55 

1485 

.007 

2243 

.010 

1300 

.006 

£ 

.60 

1328 

.006 

1892 

.009 

1047 

.005 

f- 

.65 

1119 

.006 

1659 

.009 

837 

.004 

hr 

.70 

1109 

.006 

1537 

.009 

968 

.005 

p 

.75 

1100 

.006 

1396 

.008 

725 

.004 

1 

.80 

933 

.006 

1240 

.008 

668 

.004 

.85 

919 

.006 

1202 

.008 

541 

.004 

|- 

r 

.90 

949 

.007 

1098 

.008 

494 

.004 

% 

% 

.95 

978 

.007 

1057 

.008 

476 

.004 

1.00 

1005 

.008 

1058 

.008 

702 

.006 

g,- 

1.05 

886 

.007 

1045 

.009 

852 

.007 

& 

1.10 

793 

.007 

907 

.008 

780 

.007 
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from  one  sample  to  the  next.  Therefore,  the  experimental  error  in 
this  determination  of  the  Kubelka-Munk  absorption  coefficient  of  atmospheric 
dust  is  estimated  to  be  on  the  order  of  25%  on  an  absolute  basis,  with 
a  precision  considerably  greater.  This  error  could  be  reduced  substantially 
by  performing  many  determinations  of  s*  and  averaging  the  results. 

DISCUSSION  OF  RESULTS 

The  measurements  made  here  give  the  Kubelka-Munk  absorption  coefficient 
for  atmospheric  dust,  in  units  of  cm-1.  In  practice  this  number  is 
found  to  be  nearly  equivalent  to  the  absorption  coefficient  as  defined 
by  the  Bouguer-lambert  law,  usually  differing  from  it  by  a  factor 
of  about  two  C6#  7,  8].  In  order  to  compare  the  results  of  the  work 
presented  here  with  that  of  previous  workers,  we  will  assume  the  two 
absorption  coefficients  to  be  equivalent,  so  that  a  corresponding 
imaginary  refractive  index,  n*,  can  be  computed  from  the  usual  expression: 


.i 

where  k  Is  the  absorption  coefficient  in  cm  *  and  X  is  the  wavelength 
in  cm.  In  Table  II  the  Kubelka-Munk  absorption  coefficient  and  corresponding 
imaginary  refractive  index  are  listed  for  different  wavelengths  for  several 
dust  samples. 

Fisher  [G  has  obtained  a  value  of  0.01  for  the  Imaginary  refractive 
index  in  the  visible  region  for  atmospheric  dust  collected  at  a  location 
in  West  Germany.  He  also  reported  a  samp  I e-to-sampl e  variation  of  about 
a  factor  of  3,  and  noted  that  there  was  little  dependence  on  wavelength 
in  the  visible  spectrum.  Volz  C2,  3,  ib]  has  reported  0.007  for  the 
imaginary  refractive  index  of  atmospheric  haze  particles  in  the  visible 
spectrum,  and  similar  values,  although  with  a  strong  wavelength  dependence, 
for  water  soluble  rainwater  residue. 

The  method  reported  here  gives  results  which  are  consistent  with 
those  of  other  workers.  An  exact  comparison  is  of  courss  difficult  be¬ 
cause  identical  samples  were  not  used  in  all  cases,  and  the  measurement 
errors  involved  are  not  easy  to  assess.  However,  it  can  be  concluded 
from  this  work  that  diffuse  reflectance  spectroscopy  can  give  a  useful 
measurement  of  the  imaginary  refractive  index  of  atmospheric  dust 
samples.  Inis  method  clearly  has  the  advantage  of  simplicity  and 
convenience,  since  it  consists  of  a  routine  total  diffuse  reflectance 
measurement.  Most  laboratory  spectrophotometers  are  capable  of  making 
such  measurements  using  accessories  designed  for  this  purpose. 
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